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Abstract 

In this work, it is shown that a non-zero vacuum energy density 
(the Casimir energy) for a scalar field appears in a continuous elastic 
solid due to the presence of a topological defect, the screw dislocation. 
An exact expression is obtained for this energy density in terms of the 
Burgers vector describing the defect, for zero and finite temperature. 



1 Introduction 



The importance of the Casimir effect |jT], ^ in nowadays physics can hardly 
be overemphasized. It hes at the core of the modern definition of vacuum 
energy of quantum fields. Ideally, these fields are placed in an unbounded, 
infinite Minkowski spacetime. In practice, however, they may undergo in- 
teractions, which can be viewed as suitable boundary conditions on them, 
or they may be defined in a background spacetime of non-trivial topology 
(presence of defects and boundaries) or geometry (curvature and/or torsion). 
This causes variations in the zero-point energy of the quantum fields, and 
the difference between the modified energy and the Minkowskian one is the 
so-called Casimir energy (see, e.g., Ref. for examples.). 

The Casimir effect in the presence of conical topologies has been exten- 
sively studied in the context of cosmic strings 0, |^, |], 0, but the applications 
in condensed matter systems have remained quite unexplored. Topological 
defects in solids (disclinations and dislocations) offer a natural setting for 
vacuum polarization to occur, as pointed out in a recent paper by one of us 
(F.M.) 1^. It was studied therein the Casimir effect in the presence of a discli- 
nation. The natural step further is to extend such results to the dislocation 
case, as these are more common in a solid |^. 

In this work we calculate the Casimir energy for a scalar field in the 
presence of a screw dislocation, which is done by studying this field in a 
spacetime endowed with a non-Euclidean metric simulating the defect 



The reason as to why we deal with this defect is that it proved to be the 
most tractable one after the disclination case, and also because it is related 
to the still unsolved problem of the Casimir effect in the presence of a rotat- 
ing cosmic string |Tl|, which brings naturally cosmology and sohd state 
physics together formally, for mutual benefit. The scalar field is the simplest 
field and can be immediately generalized to the electromagnetic and phonon 
fields cases by just decomposing these along their polarization axis. Also, in 
Ref. 1^ it was argued that fermionic fields (electrons in this case) should 
have a much smaller associated Casimir effect than the electromagnetic field 
due to an exponential damping related to the electronic mass. 

This paper is organized as follows: in Sec. 2 we review briefly the (- 
function regularization technique as introduced by Hawking in Ref. [1^ 



We then use it in Sec. 3 to obtain an expression for the contribution of the 
vacuum energy due to the defect. In Sec. 4 we derive an asymptotic expres- 



1 



sion for large distances, useful for practical purposes. In Sec. 5 we present 
our concluding remarks. In Appendix A we give a brief definition of the heat 
kernel of an operator and its connection with the generalized (^-function, for 
the benefit of the reader unfamiliar with these concepts. The heat kernel for 
the screw dislocation is obtained in Appendix B. Throughout the paper we 
use units in which c = h = ks = 1- 



2 (^-function regularization 



The (^-function technique [0, |13[ is a powerful tool to calculate the Casimir 
energy. Among its advantages we have that the effects due to temperature 
may be implemented in a very straightforward manner, via the Euclidean 
time formalism [|14|. As we have a spatial defect, we shall be concerned with 



an ultrastatic spacetime manifold A^, i.e., A4 = RxH, where E is modified 
by the defect. We may perform a Wick rotation of the time coordinate x^to 
the imaginary time r = ix^, and impose that the scalar field if be periodic in 
imaginary time with period P = -f , the inverse temperature. In this context, 
the partition function for the scalar field with mass m at temperature T is 
defined as (our notation and outlook follow p^ ): 

Zp= I V<^exp(-l f vpL4¥?c?^a;) = exp[^C^(0 I ^)], (1) 

where /i is an arbitrary renormalization parameter from the path integral 
measure, L4 = — + w? , and we represented by C^(0 | ^) the derivative 
at zero of the (generalized) (^-function of the operator L4 (see Appendix A). 
After a Wick rotation the topology of M. becomes x E, and A^^stands 
for the four dimensional Laplace-Beltrami operator for the (Euclidianized) 
metric on M.: 

A?! = \^9rQ.\ , (2) 

where g^J^y is the metric tensor on 7W = S*^ x S, and g is its determinant. 
For Wick-rotated ultrastatic manifolds, one has that the LB operator can be 
decomposed as A^^^ = —dr — ^lb^ now A^^^ is the LB operator for the 
metric on S (spatial part of the metric on A^). 

Eq. (1) provides all relevant regularized physical quantities via the usual 
formulas of thermodynamics, with one difference from the ordinary partition 
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function: the presence of a vacuum term, which gives the zero-temperature 
contribution. This can be isolated in the hmit (3 oo. By denoting L3 = 
—^iB + can use the relation between the heat kernel of an operator 

and its ^-function (Appendix A) and some algebraic manipulations to obtain 

H ■ 



V4vrr(s) V7i"r^"^ 



n=l 



00 

.s-1/2 



2 V47r ~, Jo 



^-(n^pf/^trj.^ e-'^^/t^'dt, (3) 
from which it is easily deduced [15,16], that 

wo, r X , 

n=l ' 

in the case one has scale independence of the Casimir energy Q In Eqs. (3) 
and (4), Tr e~*^^ stands for the trace of the heat kernel of L3 (Appendix A). 
We shall be interested in the free energy of the field F = — lnZ^//5, which 
immediately yields 

F = -C(-l/2 I L3) - ^ V / t-3/2e-M)V4t2^^ e-tL3^^_ (5) 
2 V47r „=i 



This last equation shows clearly the splitting between the T = term and 
the contribution from thermal excitations, the latter being identified with 
the ordinary thermodynamical free energy. One has hereby a conceptually 
simple method to calculate the Casimir energy, just by obtaining the heat 
kernel of L3 and its ^-function. We proceed to do this in next section, for 
the explicit case of a screw dislocation. 



3 Casimir energy around a screw dislocation 

In this section we basically use a variant of Eq. (5) to obtain an expression 
for the Casimir energy. We follow the theory of defects/three-dimensional 

"'^The Casimir energy raay depend on the scale /i and in this case the expression is 
more comphcated. In Ref. it is displayed the general formula. In our case, and in the 
majority of the physically relevant systems, one has scale independence, and therefore we 
omit further comments on this. For a discussion of these issues, see Ref. p3]. 
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gravity of Katanaev-Volovich [jT^ to model a screw dislocation by endowing 
the elastic infinite continuum representing the solid with a metric given, in 
cylindrical coordinates, by [ pi?| : 



ds'^ = dr'^ + r'^dcfP' + {dz + nd(t)f , (6) 

where < < 27r. This metric represents an infinite defect along the z-axis, 
resulting from "cutting " the space along the semiplane </> = 0, r > 0, and 
"sliding " it along the z axis by an amount 6, the modulus of the Burgers 
vector. This process gives a helix-like topology to the space, the step of 
which is K = In the case 6 = we obviously recover the flat Euclidean 
metric, describing the medium in the absence of dislocations. The metric in 
Eq. (6) is locally flat, as can be easily seen by performing the coordinate 
transformation Z' = z + ncj): 

ds" = dr^ + r2rf02 + dz'\ (7) 

which is fiat everywhere except at r = 0, where it has a 5-function singularity 
in the torsion [|lO|.One may infer the helix- like behavior by noting that the 
periodicity of induces the identification: 

it, r, 0, Z') ~ (t , r, + 27r, + 2t:k) . (8) 

As we mentioned before, our interest is in deriving the heat kernel of the 
operator L3 associated with this metric. This is done in Appendix B, and 
the result is: 



2, (r^ + r ^) 
^ — ?Ti TT 00 pcyo 

e e 4t ^ -m{0-</,') / iu{z-z') 



X 



Crr' \ 
'2t I 

where Ip{z) is the modified Bessel function of the first kind. If we consider 
r = r', Eq.(9) becomes: 

2 
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The so-called local C-function, which plugged into Eq. (5) gives the free 
energy density (see Appendix A), is obtained by setting r = r ' in Eq. (9) 
and applying a Mellin transform 



1 f°° 

C(s, r) = — - / f''K{T, f, t)dt. (11) 
i {s) Jo 

Before performing the Mellin transformation, however, we shall work Eq. 
(10) out to reshape it into a more convenient form. Let us consider the well 



known integral representation [IS]: 



1 sin UTT 

(z) = — e"=°'^ cosue dO - / e-""°'^^-''^(ix, (12) 

2ll J -TV TT Jo 

valid for Re z > Q and Re u > 0. We substitute it into Eq. (10) to obtain 

2, r2 

K(f,f,t) = , \, y [— dee^'°'' due-^ 'x 

1 f°° 2 

COS [{n + Ku) 9] / du e~'^ * sin \ n + \ x 

TT J-oo 

/ e-^™'^^e-l"+''^l^]. (13) 



We shall deal separately with the integrals. We have that [l^ : 



/oo 2 /tT K^e"^ 

dv e"*" * cos [{n + kz/) 61] = J-e ~ cos{n9). (14) 
-oo V t 

We use here the Poisson summation formula 0: 

oo oo 

F{n) = 2n ^ c(27m), (15) 

n=— oo n=— oo 

where 

1 r°° 

c{a) = — dx e'"^F (x) . (16) 
27r j-oo 

This easily gives the well-known formula: 

oo oo 

^ S(x-n)= e''""", (17) 
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which provides the following representation: 



oo 

^ COS (n6 

n=—oo 



(18) 



We recall that 6 (Xx) ~ A ^S{x), and upon substitution of Eq.(14) into 
Eq.(13), and then substituting Eq. (18) into Eq. (13), we obtain: 



K{f,f,t) 



E 



(47rt)3/2 SttH 

POO 2 

smn\n + Ku\ da; e-ir^°«hXg-|n+H^. (19) 
The next step is to calculate the local C-function (see Appendix A) : 



1 f'°° 

C(f,s) = — -/ f-'K{f,f,t)dt, 
1 [s) Jo 



(20) 



which by means of a simple representation of the modified Bessel function of 
the second kind, 



u/2 



(2^00') = dt exp (^-j- , 



(21) 



may be written as: 



C(^ s) 



1 



(47r)3/2r(s) 47r3r (s) i-00 Jo 

s-l 



dv / dx Fb {y, x) x 



r cosh 



i<'5_i(2rcosh ) \/i/2 + 777,2), 



(22) 



where we have defined 

00 

(i/, x)= sin TT I n + I e'l'^+'^'l^ 



(23) 



We may notice that for k = this last summation vanishes. Therefore, in 
the absence of the defect (6 = 0), we see that only the first term on the right 
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hand side of Eq.(22) remains. By substituting the Euchdean heat kernel, 
Eq.(B.7) of Appendix B, into Eq.(ll), we see that tliis term corresponds to 
the free space ^-function. The Casimir energy is renormahzed with respect 
to the flat, infinite Euclidean space, and we shall accordingly disregard the 
Euclidean term and thus obtain the defect ^-function: 



Cb{r, s) 



2 /"OO /"OO 

/ dv dx Fbiv^x) X 

r (s) j-oo Jo 



X 



V cosh 



s-1 



X,_i(2rcosh ( Ij ^/v'^ + w?). (24) 



As it stands, we have the Casimir energy density at zero temperature: it 
is simply E = ^Cb{r, s). We may get the full energy, valid for any temperature 
by considering Eq.(A.ll) of the Appendix A. It is clear that our problem is 
to calculate the integral : 



7 = / t~^'^K{f,f,t)e~—dt, 
Jo 



(25) 



for each natural number p. The procedure is entirely analogous to the former 
calculation, and it yields: 

-3/2 

1 roo roo .2„..u2^.^ 'r<^^^ 



du / dx Fjj (u, x) 

oo Jo 



cosh^ I + 



y2 _|_ ^2 



X 



X_3/2(2^ r^cosh^ (^Ij + ^Vl^2 + ^2). 



(26) 



Then, we finally have the free energy density for a massive scalar field in 
the presence of a screw dislocation at any temperature: 



1 1 /"OO /"OO 

^ = -,Ur,-l/2) + -Y.J_Ju I dxF,iu,x) 



xosh^ (f) 



y2 _|_ j^2 



.3/2(2^ r^cosh^ + M)IVi^2 + ^2) 



The massless case may also be otained by directly setting m = in this 
equation. The value Cb{f^, —1/2) should be understood in the following way: 
one assumes that s is large enough to make the integrals in Eq.(22) converge, 
performs the integration in terms of s and afterwards one analitically con- 
tinues the resulting expression to s = —1/2. This cannot be done explicitly 
in our case due to the complexity of the expressions. Thus we have derived 
an exact expression for the Casimir energy, but because it is a complicated 
one, in the next section we shall derive an asymtotic formula for the simplest 
particular case, namely the massless field at zero temperature. 



4 Asymptotic expression 

Our problem here amounts to see what our formulas become if we take ^ 
OO. From the physical standpoint, one expects that the qualitative features 
of the energy functional are preserved at large distances, since the defect is 
axially symmetric and infinite. On the other hand, our considering the solid 
as a continuum medium is an approximation valid for large distances from 
the defect. 

Let us now consider Eq.(lO) for the heat kernel. We set m = and and 
take the following asymptotic expression for the modified Bessel function for 

\w\^oo Q: 

/|A| (w) ^ {27iwy^^^ exp - ^ (^^ - 1 /4)) • (28) 

Now, we just set w = ^ and \ = n + ku, and substitute Eq.(28) into 
Eq.(10).The result is: 

i^(f,f,t) = -4^ 5: e-T- / due (29) 

We complete the square in the last exponencial and, defining a new variable 
u = u -\ , integrate with respect to u to obtain: 

f 2 2 

,2 OO „ t- 



i^(r,r,t) = ^(l + ^)-V2 ^ e->%^^. (30) 

°^ '-^ ^ n=-oo 
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Thus, we simply expand ~ 1 + 4^3- and apply the Mellin transform, which 
yields the local zeta function for the massless field: 

1 „2s-3 '^^ 



"a + ^r-'/\n{2s), (31) 



IGtt^F (s) ■ 

where we have substituted '^^^ = C-R. (-2)5 the Riemann zeta function 

[|T^ . Now we put s = —1/2 and pick the values of the Riemann zeta function 



and gamma function in any standard table of special functions, e.g., Ref. 



[|Tq] , then we finally have the simple expression: 

c,„„„, (-i/2.f) = (1 + '^r^ - (1 + ^r^^- (32) 

We have seen that we must subtract the Euclidean contribution in order 
to obtain the renormalized Casimir energy density: 

^ren. = ^ [C{m=0) (-1/2,^) " C(m=0) (-1/2,^) U=o]- (33) 



From this expression we obtain, for instance, the lowest order term of an 

2 

expansion in powers of \: 



^ren = — + O ( . (34) 

Note that in the units we use here the energy has units of (lenght)"^ That the 
sign is positive should already be expected, since at T = the addition of the 
defect to the background increases the energy of the field, i.e., the favoured 
configuration is 6 = 0, the fiat space without defects. Only even powers of k, 
occur, since it is immaterial from an energetic viewpoint whether the "helix" 
is right-handed (6 positive) or left-handed (6 negative). 



5 Concluding Remarks 

In the present work we give an example of the importance of quantum field 
theoretical effects in solids with defects. Albeit very simplified, our model 
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for the solid (infinite, non-magnetic and with only one defect), gives already 
a non-zero Casimir energy, which can play a role in various phenomena in 
solid state physics, for instance, contributing to the specific heat of solids. 

It would be interesting to investigate an extension of this model to mag- 
netic materials. In the continuous approximation of the Ising model, for 
instance, this amounts to add a term proportional to y?^ to the argument 
of the first exponential in Eq.(l). This model would no longer yield an exact 
result but could be handled perturbatively. 

Another line of interest is to consider classical effects of the defect. Specif- 



ically, it is known |]22[ that a classical point charge in the presence of a discli- 
nation suffers a self force induced by the defect. An identical effect would be 
expected in the dislocation case, and possibly give rise to bound states. We 
are presently tackling this problem. 

Appendix A: (^-function and the heat kernel operator 

The C-function technique as used in this paper was sistematically applied 
for the first time in Ref. |jl3|, but Hawking |jl2|] shaped it into its present 



form. We briefly introduce the subject in this Appendix. For a more detailed 
account we suggest Ref. [|12| and the comprehensive (and very complete) 



book by Elizalde et al. 

One starts with the description of a quantum physical system by a Hamil- 
tonian with adequate boundary conditions and, possibly, non-trival back- 
ground metric and topology. Mathematically, this reduces to an elliptic, 
self-adjoint, second order differential operator, with correspondent boundary 
conditions. Generally, it is not possible to calculate its spectrum directly, 
but we can still define, for such a kind of operator, the (generalized) zeta 
function of H, given by: 

Cj^(s) = Tr exp (-sin ^) ,s G a (A.l) 

In the case where we have a known spectrum, {A„} ( for simplicity, we 
assume in what follows that the spectrum is discrete), the operator e"*^"^^ 
may be diagonalized, becoming an infinite "matrix", with entries e~'^^°'*'", 
and therefore the trace in Eq.(A.l) becomes: 

C^(^) = EV- (A2) 
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Of course, if one has a zero eigenvalue, it is dealt with separately. If A„ = n, 
where n is a positive integer, and the summation runs over the positive 
integers, one has just the well known Riemann zeta function |TB|, and hence 
its name. Even if the spectrum is not known, one can still obtain information 
about the (^-function by the so-called heat equation of the operator H: 

^^F{x,y,t) + HF{x,y,t) = 0, (A3) 

where x and y are spacetime points, and t is a parameter not necessarily 
equal to the coordinate time. We assume that H acts with respect to x only, 
and that we have the initial condition 

F{x,y,0) = 6{x-y). (A4) 

Of course, Eq.(A.3) reduces to the ordinary heat equation if we have H = 
—A, the ordinary Laplacian. The function F{x, y, t) (actually a distribution) 
obeying Eq.(A.3) with condition as Eq.(A.4) is called the heat kernel of the 
operator H. 

Now let {(fn (x)} be a complete set of eigenfunctions of H, orthonormal- 
ized by the standard inner product: 

<^n (x) ^Pm {x) y/gdx = 6mn, {A.5) 

where g is the determinant of the metric tensor on the spacetime, or region of 
spacetime under consideration. Of course, they have associated eingenvalues: 

In terms of these eigenvalues and eigenfunctions the heat kernel may be 
written as: 

F{x,y,t) = ^exp(-A„t)v9„(x)v9„(y). {A.7) 

n 

The trace of the heat kernel is defined as: 

Trexp{-tH) = j F{x,x,t)^dx = ^exp(-A„t). {A.8) 

Note that F{x,y,t) is the "matrix "element of the operator exp(— ti^). 
The link with the (^-function is established via a Mellin transform [[12| : 

C/^(5) = ^y^ t'~'Trexp{-tH)dt. (A9) 
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One may also define the local ^-function : 

1 /■°° 

^Hi^^-') = Y^l t'-'F{x,x,t)dt, (AlO) 
in terms of which one may adapt Eq.(5) to obtain the free energy density: 
^ = 7;Ch{x, -1/2) - E / r^/^e-("^)V4*F(x, x, t)dt. {A.ll) 

The convenience of taking local quantities is avoiding an infinite global 
(^-f unction due to an infinite spacetime. 

Appendix B: the heat kernel for the screw dislocation 

In Section III, wc displayed the heat kernel for the operator L3 = — Ax,B + m^ 
for the metric describing the screw dislocation, Eq.(6). We also showed in 
Appendix A that the problem of obtaining the heat kernel actually boils 
down to knowing the eingenvalues and (orthonormalized) eigenfunctions of 
the respective operator, and using Eq.(A.7). The eigenvalue equation can be 
written for the screw dislocation metric as: 



\d_f d\ ^_fd__ dV _SP_ 2 
r dr \ Or \d(h dz dz^ 



A complete set of orthonormal solutions is easily seen to be: 

^n,u,k (r, 0, z) = ^ J|„+«,|(A;r)e-^"'^e^^^ {B.2) 

and its complex conjugate. \ — k"^ -\- v'^ -\- m^, where A; is a positive real 
number, 1/ is a real number, and n is an integer. J|„+„,^| (/cr) is the regular 
Bessel function. Application of Eq.(A.7) gives (the measure of integration 
for k is kdk): 

K{f, f, t) = ^ En=-o. I-oo di^ e-'^'* /o~ dk ke-'''t^-in{4>-4>')eH^-z') ^ 

{B.3) 
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Using the relation |T3]: 



e-'^'^Vp (ax) Jp {(3x) xdx = ^ ^ J /, (0^ , (5.4) 
we arrive at the form: 

(5.5) 

In order to show that it is consistent, we put 6 = to show it reduces to 
the Euchdean case. Using the identity [0: 

oo 

efK)= Y: In{x)t\ (5.6) 

n=— oo 

we obtain: 

( -^^^ 

—m^f —— 

Kb=o{r, r t) = , (5.7) 

(47rt) ' 

which in fact is the Euchdean heat kernel 12011. 
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